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(54) Method and apparatus for terahertz tomographic imaging 


(57) Providing a compositional image of an object in 
real time is accomplished by illuminating the object with 
a waveform comprising pulses of electromagnetic radi- 
ation in the terahertz frequency range and measuring 
the relative time delays of pulses reflected by the object 


to determine the positions of dielectric interfaces in the 
object. According to a further embodiment, the wave- 
form includes one or more pulses of terahertz radiation, 
each pulse having a duration in the range of approxi- 
mately 10 to 10,000 femtoseconds. 
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Description 


BACKGROUND OF THE INVENTION: 


5 1 • Field of the Invention 

The invention relates generally to electromagnetic imaging, and more particularly to imaging techniques that op- 
erate in the terahertz frequency range. 

10 ■ 2. Background Art 

Time-domain measurements of electromagnetic radiation in the terahertz frequency range have been used in the 
operational environment of spectroscopy. However, the use of terahertz electromagnetic radiation to provide an image 
of an object in real time is impractical due to the relatively long acquisition time required for terahertz waveforms. This 

15 acquisition time is typically in the range of minutes. Using pulses of terahertz radiation, one presently-existing terahertz 
("T-ray") imaging technique reduces the acquisition time of a single terahertz waveform from several minutes to several 
milliseconds, while still maintaining a reasonable signal-to-noise ratio. Examples of this technique are disclosed in U. 
S. Patent No. 5,623,145 issued April 22, 1997, to Martin C. Nuss and also in EP-A-0828 143. 

However, this technique does not explicitly utilize the time-domain nature of the terahertz waveform. Consequently, 

20 the displayed images only show the transmitted or reflected power, obtained by integrating the Fourier spectrum of the 
terahertz waveform with a digital signal processor (DSP). Such an image is not a full-volume image, which is defined 
as an image that shows the locations of compositional discontinuities within the volume of an object. Images that only 
show transmitted and/or reflected power do not reveal compositional information about an object including the locations 
of discontinuities. What is needed is an improved terahertz imaging technique which provides depth resolution and 

25 compositional information about an object in real time. 

SUMMARY OF THE INVENTION: 

Providing a depth-resolved compositional image of an object in real time is accomplished by illuminating the object 
30 with a waveform comprising pulses of electromagnetic radiation in the terahertz frequency range and measuring the 
relative time delays of pulses reflected by the object to determine the positions of discontinuities in the object. According 
to a further embodiment, the waveform includes one or more pulses of terahertz radiation, each pulse having a duration 
in the range of approximately 10 to 10,000 femtoseconds. 

55 BRIEF DESCRIPTION OF THE DRAWINGS: 

FIG. 1 is a hardware block diagram and partial pictorial representation of a terahertz imaging system constructed 
according to a preferred embodiment disclosed herein. 

FIGs. 2A, 2B, 2C ; and 2D are graphs of, respectively, a terahertz waveform used to illuminate a test object; a 
40 terahertz waveform reflected by the test object: a waveform reflected by the test object to which signal processing 
techniques have been applied: and a curve from which a refractive index profile is retrieved. 

FIG. 3A shows a prior art tomographic ray image, and FIG. 3B shows a depth-resolved slice of a test object using 
the imaging techniques disclosed herein. 

45 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 


FIG. 1 is a hardware block diagram and partial pictorial representation of a terahertz imaging system constructed 
according to a preferred embodiment disclosed herein. A terahertz transmitter 103 provides a source of terahertz 
electromagnetic radiation. Terahertz transmitter 103 is controlled by means of an enable input. The state of the enable 

50 input is related to the output of the transmitter. For example, when a pulse is present at the enable input of transmitter 
103, the transmitter is turned on, providing electromagnetic energy in the terahertz frequency range for a certain time 
duration which may, but need not, be substantially equal to the duration of the pulse. By way of an illustrative example, 
when a pulse is present at the enable input, this may result in activation of terahertz transmitter 103 for a specified 
period of time (unrelated to the pulse duration), after which the transmitter 103 automatically turns off. In the case 

55 where no pulses have been applied to the enable input of transmitter 103 for a specified period of time, the transmitter 
is turned off, and terahertz electromagnetic radiation is not generated. 

The enable input of terahertz transmitter 103 is coupled to a pulse generator 101 . Pulse generator 101 generates 
optical gating pulses that have a duration in the range of approximately 10 to 10,000 femtoseconds. In this manner, 
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the pulse generator 101 controls terahertz transmitter 103 such that the transmitter generates bursts of terahertz elec- 
tromagnetic radiation where each burst is in the range of about 10 to 10,000 femtoseconds. The electromagnetic 
radiation generated by transmitter 103 is in the form of a beam of terahertz pulses. 

'The beam of terahertz pulses generated by transmitter 103 is incident upon first reflective surface 111, which 

s reflects substantially all of the beam to a second reflective surface 113. In the present example, the first reflective 
surface 111 is flat and the second reflective surface is parabolic. Substantially all of the beam reflected by the second 
reflective surface 113 is aimed at a test object 107. The second reflective 'surface 113 may be positioned such that the 
beam of terahertz pulses are incident upon the test object 107 at a nearly normal angle of incidence, and such that 
the beam comes to a focus at or near the surface of the test object 107. In an illustrative implementation of the system 

10 shown in FIG. 1 , the beam is brought to a focus substantially at the surface of test object 107 by setting the confocal 
parameter of the second reflective surface 113 to an appropriate value related to the distance between the test object 
107 and the second reflective surface 1 1 3. When the beam is brought to a focus at or near the surface of test object 
107, this results in enhanced spatial resolution corresponding to diffraction-limited performance. 

A portion of the beam of terahertz pulses that are incident upon test object 107 are reflected by the test object. 

'5 This reflected beam of terahertz pulses is then re-collimated by the second surface 113, then captured by a third 
reflective surface 109, which directs the reflected beam to a terahertz detector 105. The beam is typically reflected by 
the test object 107 at a near-normal angle of incidence. For a test object 107 having multiple reflecting surfaces within, 
the waveform returned from the test object consists of a plurality of replicas of the waveform incident thereupon, with 
each reflecting surface returning a corresponding replica of the incident waveform. If the incident waveform is a single 

20 pulse of terahertz radiation, the waveform reflected from the test object includes a plurality of pulses of various mag- 
nitudes, polarities, and temporal distortions, reflected from dielectric discontinuities within the object. 

Refer to FIGs. 2A-2D which were prepared using a conventional 3.5-inch floppy disk as test object 107 (FIG. 1 ), 
The beam of terahertz pulses is directed by reflective surfaces 111, 113 such that the beam is incident upon, and 
reflected from, substantially a single point on the floppy disk. FIG. 2A is a graph of a terahertz pulse which is incident 

25 upon the 3.5" floppy disk, FIG. 2B is a graph of a terahertz waveform reflected by the floppy disk, and FIG. 2C is a 
graph of the waveform of FIG. 2B to which signal processing techniques have been applied pursuant to a preferred 
embodiment disclosed herein. FIG. 2D is a graph of the index profile of the object, and is obtained from the waveform 
of FIG. 2C. 

The waveform 201 of FIG. 2A may be measured at the output of terahertz detector 105 (FIG. 1 ) by replacing test 

^0 object 107 (FIG. 1 ) with a mirror or other object that reflects substantially all terahertz radiation incident thereupon with 
little or no reflection loss. In such a case, waveform 201 represents an inverted-polarity replica of the terahertz pulse 
incident upon test object 107 (FIG. 1). Waveform 201 (FIG. 2A) is inverted in polarity from the actual waveform that is 
incident upon test object 107 (FIG. 1) due to the 180-degree phase shift that occurs at the aforementioned mirror or 
other metallic objects. The small oscillations which follow a main pulse 202 in waveform 201 are a result of residual 

55 water vapor in the beam path, and do not affect system operation in any significant manner. 

Waveform 205 (FIG. 2B) is a representative waveform reflected by test object 107 (FIG. 1)as measured at terahertz 
detector 105. This waveform 205 (FIG. 2B) actually consists of a series of "distorted replicas" 207, 209, 211 of incident 
waveform 201. Each distorted replica 207, 209. 211 of incident waveform 201 corresponds to a particular reflection 
from a specific compositional discontinuity of the floppy disk. One type of compositional discontinuity is a dielectric 

40 interface. Such dielectric interfaces may occur where air meets plastic, where plastic meets air or from the surfaces 
of the magnetic recording material. The polarity and magnitude of each reflection are determined by the reflection 
coefficient at each interface, and are related to the size and sign of the corresponding index step. 

In the example of FIG. 2B, the reflections resulting from the front and back plastic covers of the floppy disk, rep- 
resented by distorted replicas 207 and 211, are clearly resolved. However, the thickness of the magnetic recording 

45 material is so small that the waveforms returned from its front and back surfaces cannot be distinguished, and appear 
as a single distorted replica 209. The incident waveform 201 hardly changes its shape while traversing test object 107 
because the plastic material of the floppy disk has little absorption and little dispersion. If other types of test objects 
are used, the reflected waveforms may be significantly altered in shape. 

FIG. 2C shows a waveform 215 that represents waveform 205 to which various numerical processing steps have 

50 been applied. These steps include a filtered Fourier deconvolution of the terahertz detector 105 (FIG. 1) instrument 
response, followed by a wavelet filtering process to remove noise. This numerical processing produces a sharp spike 
at each time delay that corresponds to the position of a reflecting interface. Thus, numerical processing helps to more 
accurately determine the positions of the various interfaces. In contrast to waveform 205 of FIG. 2B, the front and back 
surfaces of the thin (—120 mm) magnetic recording material are clearly observed in the processed waveform 215 of 

55 FIG. 2C. This is consistent with the expected resolution of Lj/2, where L c = 200 microns is the coherence length of the 
THz pulse in the intervening material. In contrast, when no other reflections are nearby, the position of a reflecting 
surface can be determined with a precision of only a few microns. Note that the amplitude of waveform 205 (FIG. 2B) 
has been scaled up by a factor of four, and that the three curves have been vertically offset along the amplitude axis 
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for clarity. In FIG. 2D, the refractive index profile of the object has been extracted from the processed waveform of FIG. 
2C according to the procedure disclosed below in connection with equation (3). 

FIGs. 3A and 3B are images of a test object obtained using the system of FIG. 1 . The test object is a floppy disk, 
and the image of FIG. 3A represents a prior-art two-dimensional image of the disk obtained by measuring reflections 

5 of terahertz pulses by the disk. More specifically the image of FIG. 3A was obtained by computing the total reflected 
power using real-time processing of the reflected waveforms with a digital signal processing device (DSP), with the 
reflected power translated into a gray scale. This is roughly similar to the original T-ray imaging introduced in U.S. 
Patent No, 5,623,145, but note that reflection geometry is used in the present invention. The plastic cover of the floppy 
disk with its various features, the circular recording disk within the floppy, as well as the metallic hub in the center of 

10 the disk can be clearly identified in the image of FIG. 3A. 

In FIG. 3B, a tomographic T-ray "slice" of the floppy disk is shown at a particular vertical position (y = 20 mm), 
indicated by the dashed line in FIG. 3A. The field amplitude of the reflected T-Ray waveforms for each horizontal (x) 
position is displayed as a function of delay in this tomographic image. The absolute amplitude of the reflected waveforms 
is translated into a gray Scale, so that each single-cycle "distorted replica" of the input wave appears as a double dark 

is stripe (refer to FIG. 2B, 207, 209, 211). The vertical axis of FIG. 3B (delay) is related to depth within the floppy disk in 
the direction of propagation of the incident terahertz pulse beam. The positions of the various parts of the floppy disk 
along the propagation direction of the THz beam can be observed clearly in this tomographic picture, such as the front 
and back cover, the magnetic recording disk, and the metallic hub of the disk (as indicated). The picture also shows 
some artifacts of the technique resulting from multiple reflections between the various interfaces, such as the features 

20 observed behind the substantially opaque hub of the disk. 

Terahertz waveforms reflected by test objects may be subjected to signal processing steps in order to more readily 
ascertain the structure of layers within the test object at any given (x,y) position on the test object. The reflected wave- 
form B(t), an example of which is given in FIG. 2B, is mathematically related to the incident waveform A(t) (FIG. 2A). 
The relationship between A(t) and B(t) may be determined by convoluting the impulse response g(t) of the layered 

25 medium: 


Bj= ifg/.aAjfe. (1) 
k=0 


In the above expression, discrete-time functions are employed. Each of these discrete-time functions is defined 
35 by digitization with a time step At (e.g., B k = B(k X M)). M is the number of samples in the digitized waveforms, and 
1024 is used in the present example. The impulse response g(t) is characteristic of the test object itself and does not 
depend upon the details of the incident terahertz pulse. The coefficient g k is the amplitude of the radiation returning 
from the object in the /c* h time slot. It is determined by the reflection off of the /c* h layer as well as by transmission through 
the preceding layers, j- 1,..., k- 1, once in each direction along the propagation path of the terahertz beam. 
40 The distance d y between two adjacent layers yand y + 1 of a test object is related to the time separation of the two 

corresponding reflections Ai = t yV; - t y by: 


d,- = — — . x Dt, (2) 
J 2 x ny v ' 


where n y is the refractive index of the medium between surfaces y and y + 1. Therefore, the temporal position of each 
reflected pulse is related to the position of the reflecting interface in the test object along the direction of propagation 
of the incident terahertz beam. 

An approximation for the impulse response function g(t) is obtained by mathematically deconvoluting the measured 
50 waveforms A(t) and B(t). Given g(t), the refractive index profile of the layered structure can then be reconstructed. In 
discretized form, the y th time step corresponds to a layer of thickness d y given by equation (2), and of refractive index 
n y given by: 


n '" ri l + r.-' W 
i = 1 


4 


BNSDOCID: <EP_0864857A1J_> 


^ ^ EP 0 864 857 A1 

where the reflection coefficients r^ are defined iteratively in terms of g^ and the previous reflections, by the relation: 


Vt 1 - 

ru = gk- II, . (4) 


J = 1 


Here r, = g^ and, by definition, the initial reference plane has n 0 = 1. These expressions are valid in the context 
10 of small reflection coefficients, when multiple reflections are negligible. Also, note that absorption and dispersion effects 
are neglected, 

A full three-dimensional representation of the object can be constructed just as rapidly as a two-dimensional image, 
because all of the information about the third dimension (the depth) is contained in a single waveform (i.e. , the waveform 
of FIG. 2B and/or FIG. 2C). For example, the tomographic image in FIG. 3B consists of 217 waveforms and took roughly 

'5 10 seconds to acquire, with each waveform recorded in 50 ms. 

Although tomographic T-ray imaging may be applied to various operational environments in which ultrasonic im- 
aging has been used, a commonly encountered problem in ultrasonic imaging is the large difference in acoustic im- 
pedance between air and liquid or solid objects, which necessitates some form of index matching. For terahertz waves, 
the dielectric constants of many dielectrics are not too different from that of air, and index matching is not required for 

20 T-ray tomography. Many materials such as plastic, cardboard, wood and rubber have good transparency in the terahertz 
frequency range. Hence, the techniques disclosed herein can be effectively used in many quality control applications, 
for example to detect voids or cracks. Another application area for T-ray tomography is burn diagnostics. While the 
strong water absorption precludes the use of THz radiation in biomedical research inside the body, it may be possible 
to obtain quantitative and highly sensitive measurements of burn depth and burned tissue properties using THz tom- 

25 ography in dermatological medicine. 

In conclusion, disclosed herein are improved terahertz imaging techniques that operate in a reflection geometry 
In the reflection geometry, timing information is correlated to depth information, and tomographic slices of objects can 
be obtained to provide 3-dimensional T-ray tomography on a wide range of test objects. In T-ray tomographic imaging, 
one can isolate successive reflections in the time domain, and extract information about each layer of a layered test 

30 object individually. In this fashion, detailed spectroscopic information can be obtained about sub-surface layers, with 
the potential for material identification. The range of potential applications for T-ray tomography is extremely broad, 
encompassing such wide-ranging fields as biomedical imaging, package inspection, and quality control. 

* - tj 

35 Claims , < 

1. A method of providing a compositional image of an object in real time including the steps of: *. : , ■- ^ 

(a) illuminating the object with a waveform comprising pulses of electromagnetic radiation in the terahertz * 
•to frequency range, and 

(b) measuring the relative time delays of pulses reflected by the object to determine the positions of disconti- 
nuities in the object. 

2. The method of claim 1 wherein the discontinuities are dielectric interfaces. 

45 

3. The method of claim 1 or claim 2 wherein the step of measuring further includes the step of using the relative time 
delays to provide a full-volume image of the object in real time, wherein the full-volume image is defined as an 
image that shows the locations of compositional discontinuities within the volume of the object. 

so 4. a method of providing a compositional image of an object in real time including the steps of: 

(a) illuminating the object with a waveform comprising pulses of electromagnetic radiation in the terahertz 
frequency range, the object reflecting at least a portion of the electromagnetic radiation as reflected pulses; 
the electromagnetic radiation traveling substantially along a direction of propagation defined by a set of parallel 

55 axes, 

(b) receiving the reflected pulses at a detector, 

(c) measuring the arrival times of each of a plurality of reflected pulses as received at the detector; and 

(d) from the arrival times, determining the location of dielectric interfaces in the object along the direction of 


BNSDOCID: <EP 0864857 A1 J_> 



10 


EP0 864 857 A1 

propagation. 

5. The method of claim 4 wherein the pulses of electromagnetic radiation include a first pulse separated in time from 
a second pulse, and the step of measuring further comprises the step of determining elapsed time between receipt 
of the first pulse as reflected by the object and receipt of the second pulse as reflected by the object. 

6. The method of claim 5 wherein the step of measuring further includes the step of using the elapsed time to provide 
a full-volume image of the object in real time, wherein the full-volume image is defined as an image that shows 
the locations of compositional discontinuities within the volume of the object. 

7. The method of claim 4 wherein the pulses reflected by the object are analyzed in the time domain to provide a 
3-dimensional tomographic image of the object. 

8. The method of claim 4 wherein the analysis of the reflected waveform uses the arrival times of reflected pulses to 
is determine the location of dielectric interfaces along the direction of propagation of the terahertz radiation. 

9. The method of any of the preceding claims wherein the pulses of electromagnetic radiation include one or more 
pulses having a duration in the range of approximately 10 to 10 ; 000 femtoseconds. 

20 10. An apparatus for providing a compositional image of an object in real time including: 

(a) a source of electromagnetic energy for illuminating the object with pulses of radiation in the terahertz 
frequency range, and 

(b) a measuring device for determining the positions of discontinuities in the object by measuring the relative 
25 time delays of pulses reflected by the object. 

1 1 . The apparatus of claim 1 0 further including an imaging device coupled to the measuring device for using the relative 
time delays to provide a full-volume image of the object in real time, the full-volume image including compositional 
information about the object. 


30 


12. An apparatus for providing a compositional image of an object in real time including: 


(a) a source of electromagnetic energy for illuminating the object with pulses of electromagnetic radiation in 
the terahertz frequency range the object reflecting at least a portion of the electromagnetic radiation as reflected 

35 pulses, the electromagnetic radiation traveling substantially along a direction of propagation defined by a set 

of parallel axes, 

(b) a detector for receiving the reflected pulses, 

(c) a measuring device coupled to the detector for measuring the arrival times of each of a plurality of reflected 
pulses as received at the detector: and 

<*o (d) a calculating device for determining, from the arrival times, the location of dielectric interfaces in the object 

along the direction of propagation. 

13. The apparatus of claim 12 wherein the pulses of electromagnetic radiation include a first pulse separated in time 
from a second pulse, and the measuring device further comprises a device for determining elapsed time between 

■*$ receipt of the first pulse as reflected by the object and receipt of the second pulse as reflected by the object. 

1 4. The apparatus of claim 1 3 further including an imaging device, coupled to the measuring device, adapted for using 
the elapsed time to provide a full-volume image of the object in real time, wherein the full-volume image is defined 
as an image that shows the locations of compositional discontinuities within the volume of the object. 

50 

15. The apparatus of claim 14 for use with an object comprising a gas or a gas mixture, the imaging device being 
adapted to generate an image of the gas or gas mixture using linear predictive coding. 

16. The apparatus of any of claims 1 2 to 1 5 wherein the measuring device analyses pulses reflected by the object in 
55 the time domain, so as to permit the imaging device to provide a 3-dimensional tomographic image of the object. 

1 7. The apparatus of claim 1 2 wherein the analysis of the reflected waveform uses the arrival times of reflected pulses 
to determine the location of dielectric interfaces along the direction of propagation of the terahertz radiation. 
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1 8. The apparatus of any of claims 1 0 to 1 7 wherein the pulses of electromagnetic radiation include one or more pulses 
having a duration in the range of approximately 10 to 10,000 femtoseconds. 

19! A method of providing an image of an object in real time including the steps of: 

5 

(a) illuminating the object with a waveform comprising pulses of electromagnetic radiation in the terahertz 
frequency range, the object reflecting at least a portion of the electromagnetic radiation as reflected pulses, 
the electromagnetic radiation traveling substantially along a direction of propagation defined by a set of parallel 

axes, 

10 (b) receiving the reflected pulses at a detector, 

(c) measuring the arrival times of each of a plurality of reflected pulses as received at the detector: 

(d) signal processing and filtering the plurality of reflected pulses to obtain an approximation of an impulse 
response function for the object: 

(e) converting the measured arrival times into depth information for the object: and 

*5 (f) obtaining one or more materials properties of the object from the impulse response function for the object. 

20. The method of claim 19 wherein the signal processing step includes deconvolving the reflected pulses with a 
reference waveform. 


20 21. The method of claim 19 or claim 20 wherein the filtering step includes a wavelet filtering process for removing 
noise accumulated in the deconvolution process. 

22. The method of any of claims 19 to 21 wherein the materials properties include a refractive index profile for the object. 

25 23. The method of claim 22 wherein information about the refractive index profile of the object is determined from the 
magnitude of the reflection coefficient of each dielectric discontinuity. 

<" r 

24. The method of any of claims 1 9 to 23 further including the step of calculating signal attenuation of reflected pulses * 5? 

due to one or more reflections and/or absorptions from compositional layers situated between the illuminating 
30 electromagnetic radiation and a specified discontinuity in the object. 

Zr ■;- -S3. 

35 ; M 
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